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Abstract — In this paper we present Ethane, a parallel search 
algorithm specifically designed for its execution on heterogeneous 
hardware environments. With Ethane we propose an algorithm 
inspired in the structure of the chemical compound of the same 
name, implementing a heterogeneous island model based in the 
structure of its chemical bonds. We also propose a schema 
for describing a family of parallel heterogeneous metaheuristics 
inspired by the structure of hydrocarbons in Nature, HydroCM 
(Hydrocarbon inspired Metaheuristics), establishing a resem- 
blance between atoms and computers, and between chemical 
bonds and communication links. Our goal is to gracefully match 
computers of different power to algorithms of different behavior 
(GA and SA in this study), all them collaborating to solve the 
same problem. The analysis will show that Ethane, though simple, 
can solve search problems in a faster and more robust way than 
well-known panmitic and distributed algorithms very popular in 
the literature. 

I. Introduction 

Metaheuristics are an important branch of research since 
they provide a fast an efficient way for solving problems. In 
many cases, parallelism is necessary, not only to reduce the 
computation time, but to enhance the quality of the solutions 
obtained. Many parallel models exist, both for local search 
methods (LSMs) and evolutionary algorithms (EAs), and even 
parallel heterogeneous models combining both methods are 
present in the literature |2J |3|. 

In a modern lab, it is very common the coexistence of many 
different hardware architectures. It has been proven that such 
heterogeneous resources can also be used efficiently to solve 
optimization problems with standard parallel algorithms Q 
||9l 1 10], but there exist few works about the design of specific 
parallel models for an heterogeneous environment. 

In this paper we propose a parallel search algorithm de- 
signed for its execution in a heterogeneous platform. We 
will also present a draft of a general model for describing a 
family of heterogeneous metaheuristics specifically designed 
for its execution in heterogeneous hardware environments, 
being inspired in the structure of the hydrocarbons that can 
be found in Nature. 

Our contribution is not only methodological, but we also 
have carried out an analysis in order to study the behavior 
of our proposal. For our analysis, we have implemented two 
versions of the algorithm making use of two well-known 



metaheuristics: steady state Genetic Algorithm (ssGA) and 
Simulated Annealing(SA). We have compared our proposal 
with the panmictic versions of these algorithms and with a 
unidirectional ring of ssGA islands executed on the same 
hardware infrastructure. Our results show that the running 
times of our proposal are faster in some cases and more robust 
in the rest than the reference ssGA ring. 

This paper is organized as follows. The next section (Section 
PI provides a brief review of decentralized and parallel meta- 
heuristics. The Section III explains the proposed algorithm 
and the model that arises from its chemical inspiration. In 



Section IV we describe some common performance measures 
which are used in this work. Section [V] contains the problems, 
parameters and infrastructure used in our study. The analysis 
of the tests is exposed in Section [Vl] Eventually, concluding 
remarks and future research lines are shown in Section I VII I 



II. Decentralized and Heterogeneous Parallel 
Metaheuristics 

In this section we include a quick review on the existing 
implementations of decentralized and parallel metaheuristics, 
as well as on heterogeneity. We also include a description of 
the metaheuristics used in our heterogeneous algorithm. 

Many parallel implementations exist for different groups of 
metaheuristics. We will focus in two of the more common 
families of metaheristics: Evolutionary Algorithms (EAs) and 
Local Search Metaheuristics (LSMs). On the one hand, EAs 
are population based methods, where a random population 
is instantiated and enhanced throw a Nature-like evolution 
process. On the other hand, only one candidate solution is 
used in LSMs, and it is enhanced by moving through its 
neighborhood replacing the candidate solution with another 
one, usually one with a better fitness value. EAs commonly 
provide a good exploration of the search space, so they are also 
called exploration-oriented methods. On the contrary, LSMs 
allow to find a local optima solution and subsequently they are 
called exploitation-oriented methods. Many different parallel 
models have been proposed for each method, and here we 
present the more representative ones. 




Fig. 1: A panmictic EA (a), and two structured EAs: distributed (b) and cellular (c). 



A. Parallel EA Models 

A panmictic EA applies its stochastic operators over a single 
population, which makes them easily parallelizable. A first 
strategy for its parallelization is the use of a master-slave 
approach where evaluations are performed in parallel but the 
population, unless divided, is treated as a whole, maintaining 
its panmictic behavior. It could be interesting for many tasks, 
but it does not offer the benefits of a structured population. 
Therefore, we are going to focus in structured populations, 
which leads to a distinction: we can distinguish between 
cellular and distributed EAs (Figure [T]i. 

• Distributed EAs (dEA): In the case of distributed EAs, 
the population is divided into a number of islands that 
run an isolated instance of the EA (Figure [TJ?). Although 
there is not a single population the sub-populations are 
not completely isolated: some individuals are sent from 
one population to another following a migration scheme. 
With this model there only exists a few sub-algorithms 
and they are loosely coupled. 

• Cellular EAs (cEA): In the cellular model, there exists 
only one population which is structured into neighbor- 
hoods, so that an individual can only interact with the 
individuals inside its neighborhood (Figure [T]:). Different 
neighborhood structures can lead to a different behavior 
With the cellular model there exists a large number of 
sub-algorithms and they are tightly coupled. 

B. Parallel LSM Models 

Many different parallel models have been proposed for 
LSMs, but there exist three models that are widely extended in 
the literature: parallel multistart model, parallel moves model, 
and move acceleration model (Figure |2|i. 

• Parallel multistart model: In this model, several indepen- 
dent instances of the LSM are launched simultaneously 
(Figure [2^). They can exchange individuals following a 
migration scheme. This model can usually compute better 
and more robust solutions than the panmictic version. 

• Parallel moves model: This model is a kind of master- 
slave model where the master runs a sequential LSM but. 



at the beginning of each iteration, the current solution is 
distributed among all the slaves (Figure |2j5). The slaves 
perform a move and return the candidate solution to the 
master, which selects one of them. This model does not 
alter the behavior of the algorithm. 

• Move acceleration model: The quality of each candidate 
solution is evaluated in a parallel centralized way (Figure 
[2};). It is useful when the evaluation function can be itself 
parallelized. The move acceleration model does not alter 
the behavior of the algorithm. 

In both EAs and LSMs parallel models, each sub-algorithm 
includes a phase for communication with a neighborhood 
located on some topology. This communication can be carried 
out in a synchronous or asynchronous manner. Many works 
have found advantages in using an asynchronous execution 
model Q. Additionally, asynchronism is essential in our study 
because of the heterogeneous hardware, which could produce 
bottlenecks, so our model communications are carried out in 
an asynchronous way. 

C. Achieving the Heterogeneity 

In the models presented above, all the sub-algorithms share 
the same search features. But we can modify the behavior 
of a parallel metaheuristic by changing the search features 
between sub-algorithms, obtaining what we call a parallel 
heterogeneous metaheuristic. Also the hardware being used to 
run the algorithm can be homogeneous or heterogeneous but 
we have not to be confused between the hardware platform 
heterogeneity and the heterogeneous software model. Parallel 
heterogeneous metaheuristics can be classified in four levels 
depending on the source of heterogeneity [ 1] : 

• Parameter level: At this level, the same algorithm is 
used in each node, but the configuration parameters are 
different in one or more of them. 

• Operator level: At operator level, heterogeneity is 
achieved by using different mechanisms for exploring the 
search space, such as different operators. 

• Solution level: Heterogeneity is obtained using a different 
encoding for the solutions in each component. 




and result 
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Fig. 2: Parallel multistart model (a), parallel moves model (b), and move acceleration model (c). 



• Algorithm level: At this level, each component can run a 
different algorithm. This level is the most widely used. 

In this paper we propose an algorithm level parallel het- 
erogeneous metaheuristic which is based in two different 
methods. We have chosen one method of each of the presented 
families, LSMs and EAs, in order to obtain a good balance 
between exploitation and exploration. The used methods are a 
Genetic Algorithm (GA) and a Simulated Annealing (SA). 

GAs are one of the more popular EAs present in the 
literature. In Algorithm [T] we can see an outline of a panmictic 
GA. A GA starts with randomly generating a initial population 
P(0), with each individual encoding a candidate solution for 
the problem and its associated fitness value. At each iteration, 
a new population P"'{t) is generated using simple stochastic 
operators, leading the population towards regions with better 
fitness values. 



Algorithm 1 Panmictic Genetic Algorithm 

Generate(P(0)); 
Evaluate(P(0)); 
t := 0; 

while not stop_condition(P(t)) do 
P'{t) := Selection(P(i)); 
P"{t) := Recombination(P'(i)); 
P"it) := Mutation(P"(t)); 
Evaluate(P"'(i)); 

P{t + 1) := Replace(P(i),P"'(t)); 
t := t+1; 
end while 



In our algorithm, we have actually used a special variant of 
the generic GA called steady state Genetic Algorithm (ssGA) 
|[T1 1 . The difference between a common generational GA and 
a ssGA is the replace policy: while in a generational GA a 
full new population replaces de old one, in a ssGA only a few 



individuals, usually one or two, are generated at each iteration 
and merged with the existing population. 

Because of its ease of use SA has become one of the most 
popular LSMs. SA is an stochastic algorithm which explores 
the search space using a hill-climbing process. A panmictic 
SA is outlined in Algorithm |2] SA starts with a randomly 
generated solution S. At each step, a new candidate solution 
S' is generated. If the fitness value of S' is better than the 
old value, S' is accepted and replaces S. As the temperature 
Tk decreases, the probability of accepting a lower quality 
solution S' decays exponentially towards zero according to 
the Boltzmann probability distribution. The temperature is 
progressively decreased following an annealing schedule. 



Algorithm 2 Panmictic Simulated Annealing 

Generate (5); 
Evaluate(S'); 
Imtiahze(To); 
k := 0; 

while not stop_condition(S') do 
S' := Generate(S',Tfc); 
if Accept(S',S",Tfe) then 

S := S"; 
end if 

Tk+i:= Update(r;,); 
k := kH-1; 
end while 



With the basis of the classic SA, many different ver- 
sions have been implemented by using a different annealing 
schedule. In our algorithm we have used the New Simulated 
Annealing (NSA) ||T2| . which uses a very fast annealing 
schedule. 




(b) 



Fig. 3: Communication schema for Ethane G (a) and S (b) 

III. Description of the Model 

In this section we present the particularities of Ethane, as 
well as we briefly outline the proposal of a generic model be- 
ing inspired in the chemical compounds called hydrocarbons. 

A. Ethane 

With Ethane we propose a Nature inspired heterogeneous 
parallel search algorithm specifically designed for it execution 
in a heterogeneous hardware platform. 

Usually, using a generic parallel model within a heteroge- 
neous platform leads to bottlenecks caused by islands with 
limited resources which can not provide good enough solu- 
tions to islands whose populations are much more evolved. 
Furthermore, many communication topologies does not take 
care of the underlying hardware architecture and could worsen 
this problem by overloading these slow islands with too much 
communication. With Ethane, we propose a communication 
schema where the most of the communication load is dis- 
tributed over the fastest nodes of the platform, and the slowest 
ones are placed as its slaves. 

The chemical compound called ethane consist of two carbon 
atoms and six hydrogen atoms, joined together with single 
chemical bonds. We have stablished a resemblance between 



® ® il 




Fig. 4: Different hydrocarbon configurations that can be found 
in Nature; their structures are the basis of HydroCM 



the atoms and the computers of the platform, and between 
the chemical bonds and communication channels. In ethane, 
each carbon atom is bonded to three hydrogen atoms, and 
there is another bond between both carbon atoms. In our 
Ethane algorithm, we propose the same schema, using two 
basic algorithms resembling different atoms, and migration 
channels resembling bonds. 

For our study, we have implemented two different versions 
of our algorithm. In Figure |3] it is shown the schema for the 
two instances of Ethane studied in this paper. Ethane G (Figure 
|3^) assigns a ssGA sub-algorithm to the central nodes, and a 
SA sub-algorithm to the slave nodes. On the contrary. Ethane 
S (Figure [3J)) allocates a S A sub-algorithm in each one of the 
central nodes, and a ssGA sub-algorithm in the slave nodes. 
With this schema, the most of the communication load falls on 
the master nodes, which are provided with the best hardware, 
moving some of the load out of the slowest nodes. 

B. An Overview of HydroCM 

From this chemical inspiration it arises a generic model 
based on the different structures of hydrocarbons. We have 
called it HydroCM (HydroCarbon inspired Metaheuristics). 
We shaped HydroCM as a generic model for a complete 
family of parallel heterogeneous metaheuristics. The goal 
of the model is to provide a schema for the islands and 
communications of the parallel algorithm to efficiently perform 
a search over a heterogeneous hardware architecture. 

Figure [4] represents some different structures for hydrocar- 
bons as we can find them in Nature. Hydrocarbons are based 
in only two different atoms, carbon and hydrogen, and each of 
them can keep a number of bounds, being one for hydrogen 
and four for carbon. 



In our model, we establish a resemblance between com- 
puters and atoms in the hydrocarbon. The bonds between 
atoms have a correspondence to communication channels, and 
double or triple bonds can be modeled as the amount of 
information being migrated or, in the case of non-population 
based algorithms, a higher migration rate. In our model, the 
fastest machines are associated with central carbon atoms 
(because of the higher computational effort caused by the 
migrations) and the slowest ones are associated with hydrogen 
atoms. 

This model provides us with plenty of different schemes for 
designing a parallel heterogeneous algorithm because of the 
amount of hydrocarbons present in Nature and their different 
architectures: linear, ring, branches... obtaining a huge amount 
of different combinations depending on the number of fast and 
slow available computers and the topology of the network. 

Ethane can be viewed as an instance of HydroCM for an 
environment composed of eight nodes, where two of them are 
more powerful than the rest, and making use of ssGA and SA. 
As well as Ethane is such an instance, we could instantiate 
many different algorithms depending on the underlying hard- 
ware architecture following the model proposed by HydroCM. 

IV. Performance Measures and Speedup 

In this section we present the performance measures that are 
going to be used for assessing the performance of the studied 
algorithms. The measures that are going to be used are the 
numerical effort, the total run time and the speedup. 

A widely accepted way of measuring the performance of a 
parallel metaheuristic is to check the number of evaluations 
of the fitness function needed to locate the optimum. This 
performance measure is called numerical effort. Numerical 
effort is widely used in the field of metaheuristics because 
it removes the effects of the implementation and the platform, 
but it could be misleading in many cases for parallel methods. 
Furthermore, the goal of the parallelism is not the reduction 
of the number of evaluations but the reduction of the running 
time. 

The most significative performance measure for a parallel 
algorithm is the total run time needed to locate a solution. In a 
non-parallel algorithm, the use of the CPU time is a common 
performance measure. While parallelizing an algorithm should 
definitely include any overhead, for example for communica- 
tion, we are not able to use the CPU time as a performance 
measure. Since the goal of parallelism is to reduce the real 
time needed to solve the problem, for parallel algorithms it 
becomes necessary to measure the real run time (wall-clock 
time) necessary to find a solution. 

Because of the non-deterministic behavior of metaheuristics, 
it is needed to use average values for time and numerical effort. 
Although 30 runs could provide us a good estimation, we have 
executed the tests 100 times in order to obtain the statistical 
analysis. 

In our analysis we will also study the speedup. The speedup 
represents the ratio between sequential and parallel average 
execution times (-E[ri] and E[Tm] respectively). 



For the speedup to be a meaningful metric, we have to 
take care of many aspects for its analysis. Because of the 
aforementioned non-deterministic behavior of metaheuristics 
it is necessary to use average times, being these times the 
wall-clock times. The algorithms run in the single and multi- 
processor platform must be exactly the same, thus panmictic 
algorithms can not be used for the analysis. The algorithms 
have to be executed until they found the solution or a solution 
of the same quality. Since in our study we are working over 
a heterogeneous platform, our reference point is the execution 
time of the program on the fastest single processor 

V. Problems, Parameters, and Platform 

In this section we include the basic information necessary to 
reproduce the experiments that have been carried out for this 
paper First we will present the set of benchmark problems 
used for assessing the performance of our proposal. Second 
we will briefly explain the parameters used within the sub- 
algorithms and the underlying hardware and software platform. 

A. Benchmark Problems 

In order to assess the performance of our algorithms, we 
have used two problems in the analysis: the Subset Sum 
Problem (SSP) |8| and the Massively Multimodal Deceptive 
Problem (MMDP) with 6 bits 0. 

The SSP problem consists in finding a subset of values V C 
W from a set of integers W = {wi, u'2, such that the 
subset sum approaches a constant C without exceeding it. We 
have chosen an instance with 2048 random integer numbers 
in the range [CIC*] following a Gaussian distribution. 

MMDP is one of so called deceptive problems. Deceptive 
problems are specifically designed to make the algorithm con- 
verge to wrong regions of the search space, decorrelating the 
relationship between the fitness of a string and its genotype. 
In MMDP a binary string encodes k 6-bit sub-problems which 
contribute with a partial fitness depending on its number of 1 's 
(unitation) following Table [l] We have used an instance with 
strings of 150 bits so that the global optimum is /c = 25. 

B. Parameters of the Algorithms and Platform 

The parameters used in every ssGA sub-population are a 
size of 64 individuals, a crossover probability of 0.8 and 
a mutation probability of 4.0 divided by the chromosome 
length. For the SA, we used the same mutation probability. 
For the SSP the chromosome length is 2048 an in the case 
of MMDP its length is 150 for both algorithms. In the case 
of the panmictic ssGA, the population size has been set to 64 
individuals because larger populations have performed much 
worse than smaller ones for the proposed problems, and they 
have been not able to find the solution of the benchmark 
problems in a reasonable time. 

We have chosen a migration frequency of 50 iterations 
for all the configurations after several initial preliminary 



TABLE I: Bipolar deception (6 bits) sub-function value 



TABLE II: Number of evaluations for the tested models and 
panmictic algorithms 



#ONES 


sub-function value 





1.000000 


1 


0.000000 


2 


0.360384 


3 


0.640576 


4 


0.360384 


5 


0.000000 


6 


1.000000 




2 3 4 

Number of ones 



experiments. The number of individuals migrated are 1 in all 
cases. For the ssGA, the emigrant is randomly selected and 
the immigrant replaces the worst individual of the population. 
In the SA, the immigrant is treated as a new move. 

The hardware infrastructure used in our analysis consists of 
8 different machines: 2 of them have an Intel Core 2 Quad 
Q9400 @ 2.66GHz processor and 4GB of RAM (namely Type 
A, fast), the other 6 computers have an Intel Pentium 4 @ 
2.4GHz processor and 1GB of RAM (namely Type B, slow). 
All the computers are managed by a GNU/Linux distribution, 
being Debian 5.0 for Type A, and SuSE 8.1, Debian 3.1 
and Ubuntu 6.10 for Type B. The computers are connected 
by a Gigabit Ethernet Network. The algorithms have been 
implemented in Java in order to support both hardware and 
software heterogeneity. For the purpose of the analysis the 
version 1.6.0_01 of the Java Virtual Machine (JVM) is used 
in all the nodes. 

VI. Tests and Analysis 

In this section we analyze the behavior of Ethane, and 
compare it with the well-known ssGA unidirectional ring. 
We have analyzed the aforementioned performance measures, 
being numerical effort, total run time and speedup, as well as 
the evolution of the fitness. 

We have implemented two different algorithms based on 
Ethane. For the first one. Ethane G, we have provided the 
Type A computers with a central ssGA island, and Type 
B computers with a SA island. For the second algorithm. 
Ethane S, the fastest machines run central SA islands and the 
slowest ones run ssGA. As we mentioned above, the migration 
scheme resembles a molecule of ethane as represented in 
Figure |3] In the parallel ssGA used as reference, the islands 
have been distributed over a unidirectional ring, placing the 
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Average 


Std. Deviation 
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Std. Deviation 


Ethane G 


146418 


174433 


1572735 


919691 


Ethane S 


202815 


198696 


708231 


430353 


ssGA Ring 


214824 


239125 


786583 


805837 


Panm. ssGA 


179792 
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Panm. SA 


81737 


93627 
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TABLE III: Time 
algorithms 



ms - for the tested models and panmictic 



Algorithm 


Subset Sum 
Average Std. Deviation 


MMDP6 
Average Std. Deviation 


Ethane G 
Ethane S 


5318 

7155 


6226 

6922 


9195 4942 
3052 1546 


ssGA Ring 
Panm. ssGA 
Panm. SA 


7453 
30008 
13300 


8107 
29387 
15443 


3194 3380 

^- * 

^- * 



most powerful computers in the first and fourth place in a 
sort of MaxSumSort |6|. As we do not know the statistical 
distribution of the data, they have been statistically compared 
with Mann-Whitney U test. 

A. Numerical Effort 

In Table [ll] it is represented the numerical effort needed to 
find the optimum for each algorithm. It can be seen that our 
proposals performed better than the panmictic algorithms for 
both problems (in the case of MMPD, panmictic algorithms 
where not even able to find the optimum in a reasonable time). 
For the SSP, both Ethane versions performed numerically 
better than the reference ssGA ring, and one of the instances 
(Ethane S) performed better even for the MMDP 

From the point of view of numerical effort, all the dif- 
ferences are statistically significant according to the Mann- 
Whitney U test. Note that also the standard deviation is better 
in our algorithms, so that its behavior is more robust. We can 
see how the panmictic SA has reached the solution with less 
numerical effort because SA is a fast converging trajectory 
method, but as we will see with the analysis of the run 
time, the time needed to find a solution is worse than for the 
studied parallel models. Since the objective of our model is 
the reduction of the total execution time let us begin with the 
analysis of a more meaningful performance metric, the total 
run time. 

B. Total Run Time 



Table III shows the average execution time of each algo- 
rithm for each problem until global optimum is reached. As 
we can see, our proposals performed clearly better than the 
panmictic algorithms for both problems (remember that the 
panmictic algorithms where not able to find the optimum for 
the MMDP) as well as better than the ssGA ring does. 



As we can see in Table III Ethane G was the best per- 
forming algorithm for the SSP problem. The Mann-Whitney 



TABLE IV: Time - ms - for the tested models in a single 
processor and its speedup 



-Ethane G -" Ethane S ssGA Ring 



Algorithm 


Subset Sum 


MMDP6 


Avg. time 


Speedup 


Avg. time 


Speedup 


Ethane G 


15995 


3.00 X 


41943 


4.56X 


Ethane S 


17817 


2.49 X 


20627 


6.76 X 


ssGA Ring 


18137 


2.43X 


21227 


6.64 X 



U test gives a p-value of 0.0412 for the Ethane G compared 
to the ssGA ring, so the difference is statistically significant. 
The average time needed for Ethane G to find a solution is 
more than 30% better than for ssGA ring. 

Ethane S was the best algorithm solving the MMDP prob- 
lem, with an average time slightly better than the ssGA ring, 
but with a much lower standard deviation, as Mann- Whitney U 
test confirms with a p-value of 0.007. The standard deviation 
of ssGA is more than twice the standard deviation of Ethane S. 
This means that the two representative instances of the Ethane 
family evaluated in this work can be more efficient and more 
robust/stable than standard sequential and distributed popular 
algorithms. 

C. Speedup 



In the Table IV we can see a summary of the execution 
time of the studied algorithms within a single processor and its 
speedup with respect to the execution in the eight processors 
heterogeneous platform. As we can see, both versions of 
Ethane have obtained a better speedup than the ssGA for the 
SSP, but only Ethane S has achieved a better speedup for the 
MMDP. 



As it is shown in Table IV Ethane G has performed better 
than the reference ssGA ring even in a single processor in 
the case of SSP. Even when its performance over a single 
processor is still better, its speedup is the best of the three 
models, however, the value for the speedup is not good for 
any of the algorithms for this problem, being the value for 
Ethane G 3 x. In the case of the MMDP Ethane G has not 
achieved as good speedup as the ssGA ring. 

Ethane S still performed slightly better than the ssGA ring 
for a single processor for both problems. Even the speedup is 
better in both cases, being the best of the studied algorithms 
for the MMDP with a value of 6.76 x. In the case of MMDP 
the speedup of the three algorithms was quite good although 
lineal speedup was not reached. 

In summary, the speedup for the MMDP was quite good 
although lineal speedup was not achieved. In the case of SSP, 
Ethane G and B have not showed a very good speedup, and 
ssGA has showed even a worse speedup. This fact could 
be explained by the huge difference among the power of 
the different hardware configurations used (remember that the 
reference point for speedup is the best performing processor). 
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Fig. 5: Evolution of fitness for SSP (time - ms - vs. fitness) 
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Fig. 6; Evolution of fitness for SSP (time - ms - vs. fitness) 



In the case of SSP, the Figure shows that both Ethane 
versions clearly outperforms the ssGA ring, converging quite 
faster. We can see how Ethane G performs even better than 
Ethane S for this problem. 

For the MMDP, Ethane S performed clearly better than 
Ethane G as we can see in Figure [6] Ethane S outperformed 
the ssGA ring, but the difference is not as large as with the 
SSP 

VII. Conclusions and Future Work 



D. Evolution of the Fitness this paper we have presented a new heterogeneous 

Figures |5] and |6] are showing the median execution of each parallel search algorithm based on the structure of ethane. We 
algorithm for each problem. have also shaped a general model for designing heterogeneous 



algorithms depending on the underlying heterogeneous plat- 
form, inspired in the structures of the hydrocarbons present in 
Nature. 

We have performed a set of tests in order to assess the 
performance of our proposal, and compared it with a well- 
known state-of-the-art model, the ssGA unidirectional ring, 
and two well-known algorithms: SA and ssGA. Our tests have 
shown that our proposal can perform better in terms of time 
and numerical effort than the reference model, and Ethane is 
even able to find the solutions in a more robust/stable manner. 
Also the speedup of the proposed models is competitive with 
that of the reference model, obtaining quite good values even 
with the huge differences between the performance of the 
computers of the heterogeneous platform. 

As future work we propose to extend and deeply analyze 
the HydroCM model, as well as to assess its performance with 
different configurations and real-life applications. Our goal is 
to offer a general model for gracefully matching computers 
of different powers to run different algorithms for efficiently 
solve the same problem, in a way that an heterogeneous 
platform does not constitute a problem but, on the contrary, 
could be used as a target platform for speciaUzed new parallel 
algorithms. 
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